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Q ■ Abstract 

1-Q . A two-fold singularity is a point on a discontinuity surface of a piecewise-smooth vector 

Cd I field at which the vector field is tangent to the surface on both sides. Due to the double 

tangency, forward evolution from a two-fold is typically ambiguous. This is an especially 
serious issue for two-folds that are reached by the forward orbits of a non-zero measure set 
of initial points. The purpose of this paper is to explore the concept of perturbing the vector 

^ ' field so that forward evolution is well-defined, and characterising the perturbed dynamics 

in the limit that the size of the perturbation tends to zero. This concept is applied to a 

cn ! two-fold in two dimensions. Three forms of perturbation: hysteresis, time-delay, and noise, 

are analysed individually. In each case, the limit leads to a novel probabilistic notion of 

^T . forward evolution from the two- fold. 

o 



1 Introduction 



Systems of piecewise-smooth, discontinuous ODEs are utilised in many areas to model phenomena 
c^' involving switching events. Except in special cases, the dynamical behaviour local to a point on 

a discontinuity surface conforms to one of three scenarios: crossing, stable sliding, or repelling. 
These are illustrated in Fig. [H In the first scenario trajectories simply cross the discontinuity 
surface. In the second scenario the discontinuity surface attracts nearby trajectories. Here 
Filippov's solution [H El [3l H] may be used to define sliding motion on the discontinuity surface. 
Sliding motion has been described in models of diverse phenomena, including oscillators subject 
to dry friction O O O |8] , relay control [H [101 E]) and population dynamics [121 [ISl [IH [ISj. In 
the third scenario trajectories head away from the discontinuity surface from both sides. In this 
case forward evolution from a point on the discontinuity surface is ambiguous or non-unique and 
a differential inclusion may be used to define a set- valued solution [161 [ITJ [H] . However, from an 
applied viewpoint this is typically not an important issue because there are no initial conditions 
for which forward evolution leads to the discontinuity surface. 
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Figure 1: Schematic phase portraits of a system of autonomous, piecewise-smooth, discontinuous 
ODEs near a discontinuity surface. The top sketches illustrate the three cases for the dynamics 
local to a point on a discontinuity surface that is generic. The bottom sketches show tangencies 
which are examples of non-generic points on a discontinuity surface. 



Points on discontinuity surfaces at which the vector field is tangent to the surface on one side 
usually represent boundaries at which the nature of the discontinuity surface changes between 
one of the three generic scenarios. For example, at the visible tangency shown in Fig. [H forward 
evolution changes from sliding motion to regular motion along the tangent trajectory. Here the 
tangency is referred to as visible because the tangent trajectory is a solution of the system. In 
contrast, the tangency shown in the bottom right sketch of Fig. [T]is said to be invisible because 
the tangent trajectory is not seen. 

Points on discontinuity surfaces at which the vector field is tangent to the surface on both 
sides are known as two-folds. Such points arise generically in systems of at least three dimensions 
[T9l|20], or as co dimension-one bifurcations in planar systems |2T]. As with repelling discontinuity 
surfaces, forward evolution from two-folds is generally ambiguous. However, unlike for repelling 
discontinuity surfaces, the ambiguity of forward evolution from a two-fold is a serious issue if 
many trajectories go into the two-fold. This is the case for the two-fold depicted in Fig. |2l 
because the area of the set of all points whose forward orbits intersect the two-fold is non-zero. 

For planar systems there are several distinct generic two-folds. These may be distinguished by 
whether neither, one, or both tangent trajectories are visible, and the relative direction of these 
trajectories at the two- fold f2T\ [22] . For Fig. HI both tangent trajectories are visible and they 
point in the same direction at the two-fold. The tangent trajectories are the orbits of the left 
and right vector fields that pass through the origin. We write them as, x=''(t) = {x"^ (t) , y'^ (t)) , 
and for simplicity assume, x^(0) = (0,0). 

In three dimensions, a two-fold for which at least one tangent trajectory is visible, may exhibit 
an ambiguity similar to that of Fig. [2] in that the forward orbits of a nonzero volume of initial 
points intersect the two- fold P^ [20] . Dynamics local to generic two- folds in three dimensions 
have been systematically classified [23| [2^ [25] . In more than three dimensions, two-folds have 



similar properties [26]. Also, it has recently been shown that two-folds arise in models of simple 
circuit systems [27], and have deep connections to folded nodes of slow- fast systems [281129] . 

However, dynamical behaviour near two-folds should only be understood with the knowledge 
that as mathematical models nonsmooth equations only represent an approximation to reality. 
Discontinuities are used to model events, such as impacts, that on an appropriately fine scale are 
not instantaneous and perhaps actually smooth, although highly nonlinear. Model inaccuracies 
are crucial at sensitive regions of phase space; indeed two-folds are points of extreme sensitivity. 

To address this issue, in this paper physically motivated perturbations of a discontinuous 
ODE system are analysed with the aim of resolving the ambiguity of forward evolution at two- 
folds. In order to obtain analytical results the scenario depicted in Fig. [2] is studied because it 
is the simplest two- fold that exhibits the serious ambiguity issue of interest. Higher-dimensional 
scenarios are left for future work. 

The following three types of perturbation are considered separately. Hysteresis - a common 
feature of control systems [30] ^ is incorporated by replacing the discontinuity surface with a 
hysteretic band of width 2e. In this context, and throughout the paper, £ > is assumed to 
be small. Second, time- delay is added. Specifically, it is supposed that the functional form of 
the equations governing forward evolution switches at a time e after trajectories intersect the 
discontinuity surface. Small time-delay is inherent in many switching systems, such as relay 
control systems [SU] , and represents the time lag between when the switching threshold is crossed 
and a change in dynamics is enacted. Lastly, noise of amplitude e is added to the differential 
equations. Noise and uncertainty are ubiquitous in real-world systems. Additive white Gaussian 
noise is used here as this is possibly the simplest manner by which noise may be incorporated 
into the system. 

With a small hysteretic band, stable sliding motion approaching the two-fold is replaced by 
rapid switching motion. Importantly, for any e > 0, forward evolution is uniquely determined in 
a neighbourhood of the two-fold. Forward orbits pass close to the two-fold and then are directed 
away from the two-fold along a path close to either x'^(t) or x^(t). We are not interested in the 




Figure 2: Dynamics local to a generic planar two-fold for which both tangent trajectories, x"*"(t) 
and x~(t), are visible and point in the same direction at the two- fold. 



precise path that orbits take, only whether they eventually "head right" , that is follow close to 
x"''(t) after passing by the two- fold, or "head left". Given, say, an interval of initial points on the 
discontinuity surface below the two-fold, for any e > 0, some fraction of forward orbits of these 
points head right after passing by the two- fold. In the limit, £ — )■ 0, this fraction approaches a 
constant value, g^y. When the system is instead perturbed by the inclusion of time-delay, the 
dynamics fit a similar description, but in the limit, e — > 0, the fraction of forward orbits that 
head right approaches a different value, ^td- The addition of noise yields a system of stochastic 
differential equations that have a unique probabilistic solution for any e > 0. From a suitable 
initial point below the two-fold, the stochastic system defines a probability that sample paths 
from this point eventually head right. As e — )■ 0, this probability approaches a value, Qno- Thus 
in each case, with e > forward evolution is well-defined, and upon passing near the two-fold 
most orbits (or sample paths in the case of the noise) follow a path close to x+(t) or x~(t). In 
each case we obtain a limiting fraction or probability, g, of trajectories that head right. 

The concept of adding noise to remove an ambiguity is not new. Under general assumptions 
the addition of noise to a system of nonsmooth ODEs with multiple solutions gives a stochastic 
system with a unique solution [3T| [32| 133] . In the zero-noise limit, this solution may represent a 
useful weighted selection of all possible solutions [311 ESj • Recently such ideas have been extended 
and applied to transport equations for fluid dynamics [36] • However, different ways of adding 
noise can give different solutions in the zero-noise limit [H7]. In addition, for a slow-fast system, 
when both the zero-noise limit and slow-fast limit are taken, the limiting solution may depend 
on the order in which the limits are applied [38] . A complete description of the zero-noise limit 
has been achieved for some simple one-dimensional systems [39] HO] . 

Variations on this concept have been applied in alternate contexts. For instance, forward 
evolution from an unstable hyperbolic equilibrium of an A^-dimensional system of ODEs is mo- 
tionless, but with arbitrarily small noise, sample paths escape the equilibrium and do so, with 
high probability, along a path near to the invariant manifold associated with the largest eigen- 
value of the equilibrium [H]. For systems of heteroclinic networks, trajectories have a "choice" 
at saddle equilibria as to which heteroclinic path to take. The addition of noise can be used to 
replace this choice with a probability [HI US] . 

The remainder of the paper is organised as follows. Equations describing the planar two-fold 
of Fig. m are introduced in ^ Perturbation by hysteresis, time-delay and noise are described 
and analysed in ^ §H and ^ respectively. Conclusions are given in ^ 



2 A planar visible-visible two-fold 

We consider the dynamics of a general two-dimensional system of ODEs near a smooth discon- 
tinuity surface. By choosing coordinates, {X,Y), such that the discontinuity surface coincides 
with the y-axis, a description of the local dynamics of the system may be written as 



X 
Y 



^W(X,F) 
/(«)(X,F) 

^(«)(x,r) 



X <0 



X >0 



(2.1) 



where we assume that the four functions that appear in (12.11) are smooth on the closure of their 
associated half-planes. In order for (12. ip to exhibit the two-fold of Fig. [2] at the origin, we must 
have 










dY 



(0,0) <0, 



^(^)(0,0)>0 



/(«)(0,0) = 

^(0'0)>0 
^(^^(0,0) > 



(2.2) 



The first two equations of (12.21) specify that the left and right vector fields are tangent to the 
discontinuity surface at the origin, the second two equations specify that both tangent trajectories 
are visible, and the third two equations specify that both tangent trajectories are directed upward. 
In order to minimise the number of parameters that must be considered, we define the scaled 
state variables 



X 



^^(0,0)^(^)(0,0) 



X, 



1 



y 



dY 



g(^){0,0) 



Y . 



(2.3) 



Since only dynamics near the origin is of interest, we expand the two halves of the vector field 
as Taylor series centred about the origin, using 0{k) to denote terms that are order k, or larger, 
in X and y. Specifically, (12. ip may be rewritten as 



X 

y 



-Ay + 0{\x\) + 0{2) 

B + 0{1) 
y + 0{\x\) + 0{2)' 
1 + 0(1) 



X <0 
X > 



(2.4) 



where 



A 



^^(0,0) 



dY 



d£ 



R) 



B 



(7(^)(0,0) 



(2.5) 



ay "(0,0) 

Consequently, A and B may take any positive values. 

The limiting fractions or probabilities, Qi^y, Qtd and g^o, determined in the the following 
sections, depend only on A and B, and not on any higher order terms. For this reason we write 
g = g{A, B) for each type of perturbation. Furthermore, in each case g satisfies the symmetry 
property 

g(j,^)=l-giA,B). (2.6) 



To see this, note that if we were to scale x and y in (12. 4p such that the vector field for x < is 
simply {x,y) = (— y, 1), plus higher order terms, then the vector field for x > would become 
{x, y) = {\y, ^), plus higher order terms. Reversing the sign of x would then recover (12. 4p . but 
with A and B replaced by -4- and ^, and the concepts of left and right switched. Therefore if 
g{A, B) represents the limiting fraction or probability of trajectories heading right, then g (-i-, -i) 
is the probability of trajectories heading left. Assuming these probabilities sum to 1 (which is 
the case for each of the three types of perturbation), we therefore have (12. 6p . 



3 Hysteresis 



In this section the switching condition of the system (12 .4^ is replaced with a hysteretic band 
spanning —e < x < e. More precisely, suppose that the equations governing forward evolution 
switch from the right half-system to the left half-system when an orbit reaches x = —e, and 
similarly the governing equations switch from the left half-system to the right half-system when 
an orbit reaches x = e. Algebraically this equates to 



X 

y 



-Ay + 0{\x\) + 0{2) 

B + 0{1) 
y + 0{\x\) + 0{2) 
1 + 0(1) 



until X = e 



until X 



(3.1) 



In [H], an analogous 0(e) hysteretic band was added to a piecewise-smooth system exhibiting a 
stable periodic orbit with sliding motion. Here the authors found that the periodic orbit persisted 
only for relatively small values of e. More generally, hysteresis may be the cause of complicated 
and chaotic dynamics [ISl HHl HT] . 

The stable sliding motion of (12. 4p that occurs at x = for y < 0, becomes rapid switching 
motion across —e<x<eioT the system (13.11) . see Fig. [31 Upon reaching the vicinity of the 
origin, some orbits head right, and the others head left. For each orbit, there is no ambiguity 
associated with its forward evolution. The purpose of this section is to precisely formulate, and 
then compute, the fraction of orbits that head right. 

There are two critical trajectories that separate orbits by whether they eventually head right 
or eventually head left. These trajectories have a tangential intersection with one of the hysteretic 
switching manifolds (x = ±e). The intersections occur at some points, (—e,y^) and (e,y^), as 
shown in Fig. [3l For our purposes it is sufficient to note that the tangency points satisfy 
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0(e) 
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0(e) 



(3.2) 



We follow the critical trajectories backward from their tangential intersections assuming that 
switching occurs at every intersection with x = ±e, as shown in Fig. [31 The values, y^, for 
A; > 1, are used to denote the successive intersections of the critical trajectories with a; = 
at which the trajectories are governed by the right half-system. By using odd values of k for 
the critical trajectory passing through (—e,y^), and even values of k for the critical trajectory 
passing through (e,y^), we have y^ < yk-i, for all k > 2. Now consider the forward orbit of 
a point, (0,y), with yk < y ^ yu-i for some k > 2, that initially follows the right half-system. 
Simply, if k is odd the orbit eventually heads right, and if k is even the orbit eventually heads 
left. 

If we treat the initial point, (0, y), as fixed, and decrease e to zero, the future of the forward 
orbit changes between eventually heading right and eventually heading left, and in the limit 
e — )■ its future is undefined. For this reason, we instead consider a large collection of initial 
points and study the fraction of points for which forward evolution eventually heads right. To 
make this precise, we ffist let x* > ehe an e-independent value chosen sufficiently small to ensure 
that global dynamics do not play a role. We then specify that an orbit "heads right" if the orbit 
ffist exits the region |x| < x*, at x = x*. We define the following indicator function for the 



forward orbit of any point, (0, y), that initially follows the right half-system: 

, . J 1 , if the orbit heads right 
^ 10, otherwise 



(3.3) 



(The assumption that orbits initially follow the right half-system is taken without loss of gen- 
erality and the value of ^hy is not affected by this choice.) For any interval of y-values, /, we 
let 

qi = ^' r : , (3-4) 



// dy 



represent the fraction of the interval / for which forward evolution eventually heads right. Already 

we have shown 

Vk^i - Vk 

q[yk,yk~2\ 



yk^2 - yk 



for /c = 3, 5, 7, . . . . 



(3.5) 




Figure 3: A sketch of phase space near the origin of the hysteretic system (13. ip . The thick 
curves are the critical switching trajectories that divide orbits that eventually head right and 
orbits that eventually head left. For instance, as indicated, the forward orbit of any point, (0, y), 
with ys < y < y2, that initially follows the right half-system, eventually heads right, whereas if 
y2 < y ^ Hi, then the orbit eventually heads left. 



Our goal is to derive 

Qhy = \iraq[y^^^,y^^^], (3.6) 

where [ymm, 2/max] is a small e- independent interval of the y-a.xis below the origin. In the following 
paragraphs we will see that g-ty is independent of ymm and i/max, and thus g-ty is an fundamental 
quantity of the unperturbed system. Before we consider such e-independent intervals, it is useful 
to first compute an explicit formula for q[y,,,y^_2]- 

For any initial point {xo,yo), with xq = 0{e) and yo = 0{y/e), if evolution (forward or 
backward) is governed purely by the left half-system up to a time t = O (y/e), then the location 
of the orbit at this time is given by 

x(^)(t) = xo - Ayot - -^e + O (ei) , 
y(^){t)=y, + Bt + 0{e). 
Similarly, if evolution is purely governed by the right half-system, then the orbit is located at 

x(^)(t)=xo + yot + V + o(£t) , 

y^''\t)=y, + t + 0{e). 
By substituting (xo,yo) = (— f^,!/^), and {x^^\t),y^^\t)) = (0,|/i) into (13.81) . we obtain 

yi = -V2y/E + 0{e) . (3.9) 

Similar substitutions using both (13. 7p and (13. 8 p lead to 



OR 

y^ = -y^^ll + —^ + 0{e). (3.10) 



In this fashion we obtain the recursion relation 



Vk = -Jyl-2 + 4 (l + ^) ^ + ^(^) , for A; > 3 . (3.11) 

Then it is easy to show that the combination (I3.9p - (l3.1ip . gives the explicit formula 

-v^7fc+(^v^ + 0(e), forA; = l,3,5,... 

yk ^ \ I . \'-'-^^) 

-V2jA;-l + f v^ + 0(£) , for fc = 2,4,6,... 



By substituting (I3.12p into (13.50 . and considering large /c, we obtain 

A 



% 



yk,yk^2\ 



A + B 



+ O f ^ j + O(v^) , for A; = 3, 5, 7, . . . . (3.13) 



Naturally we would like to use (I3.13P to evaluate (13. 6p . Let k^in and /^max be the smallest and 
largest odd integers for which yu G [ymm,2/max]. Then, by (JS]!]), 

.[.......] = ^^^^^-"^'^-r-^'"^""'^ + 0{e) , (3.14) 




Figure 4: The fraction of orbits that head right after passing near the two-fold for the hysteretic 
system ([31]) in the hmit e ^ 0, ^hy(^, B) = ^, ( ISTT]) . 



where the error term corresponds to the 0{e) differences, ymax — l/femin ^^^ Vh 
quotient, f l3.14p . may be written as a convex combination of, q[y,,,y^_2] 



The 



max Vrn'm- 

y^-'^~y^ for all odd k 



with kr, 



2 < k < kraax- For this reason we have the bound: 



yk-2-yk ' 



(3.15) 



Z <. k "^ fcmax- -tlOWeVer, /Cmim fCmax ^ 



- ; "^min) "'max 



where the extrema are computed over all odd k with /cmi 
oo as £ — 7- 0, and since it has not been shown that the 0{^/£) term in fl3.13p remains bounded as 
/c — 7- oo, equations f l3.13p and f l3.15p are insufficient for us to derive the anticipated and intuitive 
result, f l3.16p . given below. We have overcome this technical difficulty via an involved procedure 
of carefully bounding error terms. We state the following result and provide a proof in Appendix 



Lemma 1. Given any — 1 ^ ymin < ^max < 0, 



9[S/min 



^inaxj 



-)■ 



A 



A+B 



as £ — )■ . 



(3.16) 



The condition, — 1 ^ ymin, ensures that global features of the system do not influence the 
dynamics. By (13. 6 p and (I3.16p . 



^hy = Qhj{A,B) 



A 



A + B 



(3.17) 



As mentioned above, ^hy is independent of the values of i/mm and ^max- A plot of ghy{A,B) is 
shown in Fig. HI 



4 Time-delay 

Now suppose that the switching condition of (12 .4^ involves a time-delay of length e: 



X 

y 



-Ay + 0{\x\) + 0{2) 

B + 0{1) 
y + 0(1x1) + 0(2) 
l + 0{l) 



x{t 
x{t 



£) <0 
e) >0 



(4.1) 



The inclusion of time-delay may destroy periodic orbits with sliding |44] . Dynamics and bifurca- 
tions specific to general piecewise-smooth ODE systems with time-delayed switching conditions 
are discussed in [18], and in [19] for systems with symmetry. If both hysteresis and time-delay 
are incorporated in the switching condition, there are yet more potential complications [50] . 

About the negative ?/-axis, dynamics of the time-delayed system (14. ip switches at curves, 
r~ and F"*" (shown in Fig. [5]) that may be thought of as the result of evolving the ?/-axis by 
a time e under each half-system. As in the previous section we use series solutions to obtain 
analytical results. If Xq = 0{e^) and i/q = 0(e), then evolution up to a time, t = 0(e), for the 
left half-system of (14. ip is given by 



X 



(L) 



it) 



xo - Ayot — t^ + 0{e^) 



(4.2) 



y(^)(t) = yo + Bt + 0{e') 



and for the right half-system of (14.10 . is given by 



x^ 



0(r 



/''\t) = yo + t + 0{e'). 



(4.3) 



From (I12D and (gSD we obtain 



[x,y) 



X 



ey e 

^ 2 



0{e' 



[x,y) 



X 



-eAy + 



— £2 + 0(£3) 



y 



0{e) 
-0{e, 



(4.4) 
(4.5) 



As shown in Fig. [5l upon passing close to the origin, most forward orbits cease switching 
and either head right or head left. However, there are two critical orbits that switch across 
the positive y-axis indefinitely, or at least until global features of the system induce different 
dynamical behaviour. As in the previous section, we let yk for k>\ denote the intersections of 
these orbits with the negative y-axis at which the orbits are governed by the right half-system, 
see Fig. O With the definition (13.41) . we again have 



%fc: 



yk-2\ 



yk-i - yk 
yk-2 - yk 



for k = 3,5,7, .. . 



and our goal is to calculate 



^td= Jim?[?/min,J/n.ax] , 



(4.6) 



(4.7) 
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Figure 5: A sketch of typical dynamics near the origin of the time-delayed system (14. ip . 



for some — 1 ^ ymin < 2/max < 0. The remainder of this section consists of a derivation of 
formulas involving the switching points of the two critical orbits, a recursion relation for y^, and 
a numerical evaluation of f l4.7p . 

The orbit of (14.10 that is shown in Fig. [6]is located at (0, ^o) at t = and satisfies the following 
two assumptions. First, it is assumed that for all t G (— £, 0), the orbit lies to the right of x = 0. 
Consequently, until t = e, the orbit is governed by the right half-system. Second, we assume that 
2;o < is sufficiently small that at t = e, the orbit is located at a point (xi, Zi) with xi > 0. By 
(14.30 we have 

xi = zoe + -e^ + 0{s^) , (4.8) 

zi = zo + e + 0{e^). (4.9) 

The points, {xk,Zk), for A; > 1, are used denote successive switching points of the orbit. We let 
Ti denote the time taken for the orbit to first return to x = 0, then 

Ti = -2zo + 0{e^) . (4.10) 

Similarly we let Tj,, for each k > 2, denote the time taken for the orbit to reach x = from each 
{xk-i,Zk-i). Finally we let Sk, for each k > 1, denote the elapsed time between the consecutive 
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Figure 6: A sketch of typical dynamics near the origin of the time-delayed system (14. ip . 



switching points {xk,Zk) and {xk+i,Zk+i). Note that, 5*1 - the time taken for the orbit to travel 
between {xi,zi) and (x2,Z2) - is, by definition, also equal to the length of time that the orbit 
spends in the left half-plane between (0, ^o) and (xi, zi), that is. 



Si = T^. 



(4.11) 



Similarly 5*2 is equal to the length of time that the orbit spent previously in the right, that is 
S2 = e — Ti + T2, and for k > 3, Sk = Sk^2 — ^fc-i + T^. From this, also by (14. lip , we obtain the 
simpler recursion relation: 

Sk = e- 5fc_i + Tfc , for A; > 2 . (4.12) 

Elementary use of the series solutions, (14. 2 p and (14. 3p . leads to 



n 

Zk 



h{^4-i + 'i^k-i-zk-i)+0{e') 

Xk-l - AZk^lSk^l - ^^Li + 0(£3) 
Zk-i + BSk^i + Oie^), 



for k 



2,4,6, 



(4.13) 
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Figure 7: The fraction of orbits that head right after passing near the two-fold for the time- 
delayed system (14. ip in the limit 5 — ;■ 0, gtd{A,B), computed numerically as described in the 
text. 



and 



Zk 



zl^^ - 2xk-i - Zk_i + 0{e^) , 

Xk-i + Zk-iSk-i + I^Li + 0{e'^) 
Zk-i + Sk-i + 0(e2) , 



for A; = 3, 5, 7, 



(4.14) 



For any values of A and B, the above recursion relations may be used to numerically identify 
the value, zq = yi{A,B), plotted in the inset of (|5]), for which, in the limit e — )■ 0, the orbit 
switches across the positive y-axis infinitely many times. The analogous value, yi < 0, for the 
orbit that initially heads right, then switches indefinitely, may obtained in a similar fashion. 
Alternatively, by symmetry, 

y,{A,B) = y,(^^,^'^ . (4.15) 

Note that it appears to be not possible to obtain an explicit expression for yi {A, B) from the 
recursion relations. 

For backward evolution on the negative y-axis, illustrated in Fig. [5l the following recursion 
relation may be derived from (14. 2 p and (14. 3 P (in the same manner as in the previous section): 



yk = -^\A + B 



Also 



A^(a+B+ yjyl_2/e^ + A{A + i5)) +{A + B)B ] e+0{e^) , for A; > 3 



2/2 



A 



A + B+ \A^ yl/e^ + {A + B)B\e + 0{e^) . 



(4.16) 
(4.17) 



The formulas are sufficient for us to perform an efficient numerical evaluation of the fraction, 
^[yk,yk^2]- To do this, for any A and B, yi is computed via an iterative scheme that identifies a 
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value of zq for which recursion according to fl4.12p . f l4.13p and f l4.14p . and starting from (14. 8p - 
(14. lip , gives values {xk,Zk) that lie on successively opposite sides of x = for as many values 
of k as one chooses to compute. Then (14.151) is used to repeat this process and compute yi. 
Successive i/k are computed via (I4.16P and (I4.17p . and finally (l[y^.y^._2] is evaluated by (14. 6p . 

It is expected that ^td = ^^^^k^oo^^^^e^oQ[yk,yk-2]y since this identity holds for the hysteretic 
perturbation of the previous section (it is an elementary consequence of Lemma[T]). However, since 
here an exact expression for y^ is not known, and in view of the complexity in proving Lemma 
[H a verification of this claim is left for future studies. Fig. [7] shows a plot of Qtd approximated 
by evaluating q[y^.,y^._2] at suitably large k and small e via the iterative scheme described above. 



5 Noise 

This section concerns the following formulation of (12.41) in the presence of noise: 



dx{t) 
dy{t) 



-Ay + 0i\x\) + 0i2) 

B + 0{1) 
y + 0{\x\) + 0{2) 
1 + 0(1) 



X < 



X > 



dt + e 



du 
d2i 



di2 

d22 



dWi{t) 
dW2{t) 



(5.1) 



Here e denotes the overall noise amplitude and dij are constants that control the relative strength 
of the two-dimensional noise in the x and y directions. Alternative formulations include coloured 
noise, spatially dependent noise, and multiplicative or parametric noise, each of which may be 
more suitable in certain applications. 



5.1 General aspects of the zero-noise limit 

Equation (15. ID is an example of a multi-dimensional stochastic differential equation, 

rfx(t) = 0(x) dt + sD dW{t) , 



(5.2) 



for which the drift, 0, may be discontinuous. If is measurable and bounded and the diffusion 
matrix, D, is non-singular, then (15. 2p has a unique strong solution from any initial point [311 [32| 
[511 [52]. In following discussion, Pei'x, t; xq) denotes the probabihty density function of x(t), given 
x(0) = xo. 

With e = 0, ([E2]) reduces to the ODE 

i = 0(x) . (5.3) 

If (j) is smooth, (15. 3|) has a unique solution, call it Xdet(^;xo). Straight-forward asymptotic argu- 
ments are sufficient to demonstrate weak convergence: Pe(x, t; xq) — )■ (5(x — Xdet(^; xq)) as e — )■ 0, 
where 6 denotes the Dirac-delta function [53],[5l]- If 4> is continuous but non-differentiable, (15.31) 
may have multiple solutions and p^ may converge to a weighted sum of delta functions centred 
at some of these solutions. In particular we may have 



Pe(x,t;xo) ^ (1 - ^)(5(x-Xi(t;xo)) + ^(5(x - X2(t; Xq)) 
as e — !■ 0, where xi(t; xq) and X2(t; xq) are two different solutions to (15.31) . 



(5.4) 
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In the short article |39], stochastic calculus methods are used to prove (15.41) in one dimen- 
sion when (p is continuous, odd, nondecreasing, convex on [0, oo), and with a finite value of 
limA^o /o ^^ dx, when xq = 0. Here, f? = | by symmetry, and of the uncountably many solu- 
tions to (15.31) . the Xi and X2 that appear in (15.41) are the two solutions for which x{t) ^ for all 
t > 0. These are called the minimal and maxim,al solutions because, for any other solution, x{t), 
for all t > 0, Xi{t) < x{t) < X2{t). In [ID], similar results are obtained when is not odd and the 
value of g is derived using first passage concepts. Similar results for multiple dimensions have 
recently been described in [55] . 

For discontinuous </>, it is necessary to consider generalised notions of a solution to (15. 3p [TS] . 
In [511 was is shown that solutions to (15.21) may limit only to Filippov solutions of (15. 3p as the 
noise amplitude is taken to zero. However, in the case that there are multiple solutions, this does 
not tell us which Filippov solutions contribute to the limiting density. On stable sliding regions 
there is a unique Filippov solution, and if sliding occurs for the duration of time interval, [0, t], it 
may be shown directly that Pe(x, t; xq) limits to a delta-function centred at the Filippov solution 
as e — >■ 0, via an asymptotic expansion of p [FT] . 

If is discontinuous, then, to the author's knowledge, only when is one- dimensional and 
piecewise- const ant is an exact expression for ^^(x, t;xo) known. Specifically, for any constants 
0^ G M, the transitional probability density function of 



dx{t) = { ;^+ ' ;^ < U \dt + e dW{t) , x(0) = , (5.5) 



is given by 



e^ Jo Jo ^27re2(t_n)3 V27ze^u3 ' - 

see [581 159] . Moreover, if 0~ < < 0^, then multiple Filippov solutions emanate from x = 0, 
and by using Laplace's method [60] to evaluate (15. 6p asymptotically, it may be shown that 

Pe{x, t; 0) ^ ^~^^_ 6{x - <P~t) + ^_^'^^_ 5(x - 0+t) . (5.7) 

This matches the results for continuous 0, in that cj)~t and 0"'"t are the minimal and maximal 
Filippov solutions to (15. 5p with £ = 0. Large deviations for the case (0~,0^) = (—1,1) and 
generalisations are described in [6T] . 



5.2 Asymptotics of a first passage representation 

Our interest is in forward evolution of ( 15. ip from points on the negative y-axis. For small £ > 0, 
upon passing near the two-fold, sample paths follow close to either x"'"(t) or x~(i) with high 
probability. This is illustrated in Fig. [HI Here the probability that a sample path will eventually 
head right or head left is determined via first passage methods for arbitrarily small e > 0. 

We restrict our attention to an e-independent rectangle, [— x*,x*] x [— y*,y*], in the {x,y)- 
plane, and assume x*,y* > are sufficiently small that global features of the dynamics may 
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be ignored. For any initial point, (x(0),y(0)) = {xo,yo), inside the rectangle, we let Qe{xo,yo) 
denote the probability that first passage by (15.11) to the boundary of the rectangle occurs at a 
point with a; > 0. We then define 



^no = \imQ^{0,yo;x*,y*) 

£-S>0 



(5.8) 



for, —y* < yo < 0, as this represents the probability that forward evolution from the negative y- 
axis heads right after passing by the origin in the zero- noise limit. Below it will become apparent 
that this limiting probability is independent of the values of x*, y* and yo as implied by (15. 8p . 

As described in stochastic differential equation texts, such as [52], Qe is the C^ solution to 
the elliptic boundary value problem comprised of the backward Fokker-Planck equation for (15. ip : 






-^yo + 0(|xo|) + 0(2))g 



Oi\xo\) 



0(2))^ 



d'Qe 



+ -2(di. + dt,),^. 



and the boundary conditions 

Qe{xo,±y*) = X[o,oo)(a;o) 



dxo 

e{dud2i 

Qe{ 



(1 



(^ + 0(1))^ 
f 0(1))^ 



dyo 



diod' 



12"22 



dyo ' 

dxodyo 2 



zo < 
xo > 



21 + "22. 






'X ,yo) 







Qe 



.X ,1/0 J 



(5.9) 



(5.10) 



The boundary conditions follow from the observation that first passage occurs in zero time for 
any point on the boundary of the rectangle. Therefore, as stated in (I5.10p . on the boundary, 
Qe = 1 if Xq > 0, and Q^ = otherwise. 

A description of Qe for small e may be obtained from an asymptotic expansion of (I5.9p - (l5.10p . 
Since only lim^^o Qe is of interest, only the leading order term in this expansion is required here. 
A formal justification of the expansion is a difficult analytical task deferred for future work. 




Figure 8: Six typical sample paths of (15.1 
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Here the scaling laws of the asymptotic expansion are obtained in an intuitive fashion via a 
blow-up of the dynamics of (15. ip about the origin. Specifically, with the general scaling, 



X 



X 



y 



y 



where each A, > 0, (15.11) becomes 



dx{t) 

dm 



~£-^^+^^+^»Ay{i) + O {e 



t 

-Ai+A3+min(Ai,2A2)^ 



(5.11) 



e' 



-A2+A3 



-Ai+A2H-A3„~ 



B + 0{e 



-A2+A3+mm(Ai,A2) 



) 



-A2+A3 



y(i) _)_ Q (£--Al+A3+min(Al,2A2)^ 



0{e 



-A2+A3+min(Ai,A2) 



+ 



^1-Ai + ^ 



-1-A2 + ^ 



dii 

2 ^21 £ 



£1-^1 + ^(^12 



l-A2 + 4f 
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) 

dWi{i) 
dW2{i) 



X <0 
X > 



dt 



(5.12) 



The appropriate values of Aj are those for which three terms in (I5.12p are 0(1) and all other 
terms are of higher order. This gives two possibilities. First we may have Ai = A2 = A3 = 2, but 
this is merely the scaling that in general sets drift and diffusion to the same order and thus covers 
a time-scale too short to describe stochastic departure from the discontinuity surface. Hence, we 
use the second possibihty, which is 



Ai 



4 
3 



A, 



2 
3 



A, 



2 
3 



(5.13) 



Since (I5.13P scales x to a higher degree than y, noise in y is of higher order than noise in x. By 
absorbing the magnitude of the noise in x into e, we may assume 



'^11 "•" '^12 — -'- ' 



so that the noise in x is normalised. In the limit, e — )■ 0, we define 

X— )-u, y ^ V , t— )-s, 

with which the stochastic differential equation, (I5.12p . takes the reduced form 



(5.14) 



(5.15) 





• 


" -Av ' 


du{s) 
dv{s) 


V 


V 

1 





u<0 
u>0 



ds + 



dW{s) 




(5.16) 



We now return to the boundary value problem for the probability, Q^- In terms of the 
limiting scaled variables, (I5.15p . the leading order term for the asymptotic expansion of Qe, call 
it Q, satisfies the backward Fokker-Planck equation of (15.161) : 







-^-og 



+ BIP- 

dvo 



„ dQ 



dQ 

dvo 



+ 

id'^Q 
2 9^ 



la^Q 

2 9^ 



Mo < 
Mo > 



(5.17) 



Since the boundary conditions for Q^ (I5.10p are at the boundary of a fixed rectangle, the boundary 
conditions for Q are at infinity. However, by the nature of the reduction, ( 15.17P is parabolic and 
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consequently we cannot impose a boundary condition on Q in the limit vq — ?■ — oo. This is 
not problematic because in the w-direction (15.161) has positive drift and no noise, thus for sample 
paths to (15. 161) . v{t) — )■ oo as t — )■ oo with probability 1. Q satisfies the remaining three boundary 
conditions, 

lim Q{uq,vq) = X[o,oo)(%) , lim Q{uo,Vo) = , lim Q{uo,vo) = 1 . (5.18) 

The goal is to determine ^no (15. Sp . which may now be written as g^o = ^i^n^vo^-oo Q(0, t'o)- 

5.3 The limiting probability, ^no 

The partial differential equation, (I5.17p . takes a more familiar form if we introduce the pseudo- 
time 

r = -vo . (5.19) 

In terms of r, the parabolic boundary value problem is 






(5.20) 



lim Q{uo,r) = xio,oo){uo) , lim Q{uo,r)=0, lim Q(mo,'^) = 1 • (5.21) 

r— >■— oo mq— >■— oo mq— >oo 

Reinterpreting Q as a function of uq and r, we have 

^no = QnoiA, B) = lim Q(0, r) . (5.22) 

r— >oo 

Due to the discontinuity at mq = and the presence of r on the right hand side of (I5.20p . it 
appears to be not possible to obtain an explicit expression for Q from (I5.20p - (l5.2ip . Fig. [9] shows 
a typical plot of Q{uo,r) as computed by finite difference approximations. Fig. [TO] shows a plot 
of QnoiA, B) computed by then approximating the limit ( I5.22P at a suitably large value of r. The 
following lemma provides the value of g^o at some special values of A and B: 

Lemma 2. The limiting value, i^5.22\) . of the parabolic boundary value problem, Ii5.20\) - f5.21\) . 
satisfies 

QnoiA, A') = ^^, (5.23) 

ft.o(Al) = ^, (5.24) 

\jm QnoiA, B) = ^. (5.25) 



These values are highlighted in Fig. [101 A derivation of (I5.23p -( l5.25p is given in Appendices 
iBllDl Since A and B may take any positive values, it follows from (I5.23p that all values in (0, 1) 
are possible for Qno- 
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Q{uo,r) 




Figure 9: The function Q{uo,r) when A = 2 and i3 = 4 as computed by numerically solving 
the boundary value problem f l5.20l) - fl5.2ip . For these values of A and B, by fl5.22p and fl5.23p we 
have Q{0, r) — )• ^ as r — )• oo. 








Figure 10: A plot of gno{A,B) computed by numerically solving the boundary value problem 
fl5.20p -f l5.2ip and evaluating Q{0, r) at a large value of r to approximate f l5.22p . The solid curves 
are the theoretical values of gno{A, B) as given by Lemma [21 



Finally, consider the transitional probability density function for (15. ip . call it pe(x, t;xo). 
Write xo = (xcyo), and suppose xq = and — 1 ^ yo < 0. Let tsiide denote the time taken 
for the deterministic sliding trajectory of (12.41) to travel from xq to the origin. Then for small 
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t > 4iide, the above results imply 

Pe(x,t;Xo) -^ (1 - ftio)<^(x-X"(t-tslidc)) + &o<^(x-X+(t- tslidc)) , (5.26) 

as e — )■ 0. 

6 Conclusions 

The two-fold of Fig. [2] has two important properties: (i) many orbits reach the two-fold, and, (ii) 
forward evolution from the two-fold is ambiguous. To clarify the second property, the two-fold 
has uncountably many distinct forward orbits in the sense of Filippov. Specifically there are two 
tangent trajectories, x"'"(t) and x~(t), and there are orbits that travel along the positive y-axis 
for some time, then enter and evolve in either the left or right half-plane. Local dynamics are 
described by the two-dimensional ODE system, (12. 4p . which has two parameters, A and B, that 
relate to the strength of the vector field for x < 0, relative to that for x > 0, in the x and y 
directions respectively. 

In this paper three different 0{e) perturbations to fl2.4p have been studied: hysteresis, time- 
delay and noise. Each perturbation removes the ambiguity of forward evolution and in each 
case forward evolution close to the two-fold has been analysed. This microscopic analysis of the 
two-fold uncovers a macroscopic view of the dynamics. For small e > 0, most orbits (or sample 
paths in the case of noise), follow close to either x"''(t) or x~(t). In the limit, e — )■ 0, almost 
all orbits or sample paths follow exactly one of the these trajectories. Thus the possibility of 
evolution along the positive y-axis for a nonzero length of time may be excluded. This observation 
mirrors theoretical results for stochastically perturbed continuous, non-differentiable systems in 
the zero- noise limit [391 HO] • 

With a hysteretic band of width 2e the system is described by (13.11) . with a time-delay of 
t = e the system is described by (14. ip . and with 0{e) noise the system is described by (15. ip . It 
is interesting that the three different types of perturbation exhibit three different scaling laws. 
Specifically, for hysteresis, x = 0{e), for orbits approaching the two-fold about the negative 

y-axis. For time-delay, x = O(e^), near the two-fold. Finally for noise, x = O ieaV near the 

two-fold. 

For perturbation by hysteresis and time-delay, ^hy and ^td, respectively, represent the fraction 
of orbits that head right (that is, follow close to x"'"(t) after passing by the two-fold) in the limit 
e — )■ 0. Similarly, for the stochastic perturbation, ^no represents the limiting probability that 
sample paths head right. As a function of A and B, an explicit expression was only obtained 
for f)hy, specifically (I3.17p . For gtd and g^o efficient algorithms were established to evaluate these 
values numerically. Lemma [2] gives special values of g^o- Plots of each g{A, B) are shown in 
Figs, m [7] and [TOl Each plot has the same scale and viewpoint so that they may be compared 
fairly. For example, we notice that as ^ — >■ 0, for hysteresis and time-delay ^ — )■ 0, but for noise 
^ -T- |. Also gno{A, B) appears relatively flat, but, by (I5.23p . may take value in (0, 1). 

The primary goal of this work has been to resolve the ambiguity of forward evolution from 
a two-fold. To conclude, we may interpret forward evolution of the unperturbed system (12. 4p 
from the two-fold as probabilistic, with evolution along x"'~(t) with probability, g{A,B), and 
along x~(t) with probability, 1 — g{A,B). However, the value of g depends upon the type 
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of perturbation under consideration. This observation indicates that such a regularisation of 
two-folds is potentially futile in the absence of additional information about the system. As a 
mathematical model the system represents only an approximation to real-world behaviour. For 
a given application, physical reasons may justify the use of one of the three types of perturbation 
over the other two. It remains to apply a probabilistic representation for forward evolution 
from a two-fold to a physical system. Certainly by discounting Filippov solutions that traverse 
part of the repelling sliding region, the two-piece probabilistic solution proposed here is more 
sophisticated than permitting all possible Filippov solutions. It remains to extend the ideas 
described here to two-folds in three or more dimensions that arise generically in discontinuous 
systems of ODEs. 

A Proof of Lemma [1] 

Here we write fl3.1ip and fl3.12p as 



Vk 



lyt, + 4(1 + ^6 + G{yk-2; e)e , for /c > 3 , (A.l) 



-V2Jk+^!^^ + Ck{e)e, for A; = 1,3,5, . . . 
Vk = \ ^ I — ■ lA.ij 



-v^^/A; - 1 + f yi + Ck{e)e , for A; = 2, 4, 6, . . . 
There are two main parts to the proof. First it is shown that there exists Mi G M such that 

\Ue)\<kMi, \lk<e-\ (A.3) 

given < 6* < 1. Second it is shown that 

|g[,.+4,y.+2] - 1[y.+.,yu] I < M^e'^' , VA; > e~' , (A.4) 

for some M2 G M. Lemma [1] is then a consequence of a few additional observations. The key 
point is that the error in flA.4l) is of higher order than e which enables us to usefully evaluate 

We restrict our attention to y G [y*,0], £ G [0,£:*], where y* and e* are suitably small. 
For this reason we may assume that the error in flA.ip . G, is a smooth bounded function and 
let Ml = 2max(|G|) + max(|Ci|) + |max(|C2|). Trivially ([Ql) holds for A; = 1 and A; = 2. 
Mathematical induction is required to verify (]A.3P for all k < e~^ . Combining ( jA.ip and (]A.2p 
leads to 



a+2 = G'(yfc) + --=yA; + 2+ ^ ^' 



- ^y A: + 2 + ^-^ - V2^k + ^-^C.v^ + ^ae , (A.5) 

for odd values of k, and a similar formula for even k. By appropriately bounding the large 
square-root term in (lA.Sp . which requires assuming e* is sufficiently small that 

^^^'^'^ < ^, ^ H^ , (A.6) 

2v^(l + |)Mi 
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it may be shown that 



(fc^i)g 



i4« - 41 < iG(,.)i + 141 1 1 - /^ " + ^n^+})asre ^ (^^,, 



t + 2 + tt±ia' *' 



lk+ 



(fc-i)£ 



By further using (IA.6p . noting -f > 1 ; f =, and applying the induction hypothesis, we 

obtain 

ia+2 - Ck\ < \m^ + Ml + ^Mi = 2Mi . (A.8) 

Consequently, |Cfc+2| < (^ + 2) Mi, for odd values of k. The result for even k may be shown in 
the same manner. This completes the inductive demonstration of (1A.3I) . 

Since (13 .14^ involves only odd values of /c, for remainder of the proof it is assumed that k is 
odd (even values of k may be treated similarly). Substituting (IA.2P into (13.51) yields 



^+o(i) + ^y^i + 5(a_i-a) 

<i[yk,yk-2\ = — I • (A.9) 



By applying (1A.3I) . an evaluation of (1A.9I) at A; = 0(£: ^) gives 

%....-.] = 7^ + 0{e') + O (.^(^-^^)) . (A.IO) 

To minimise the magnitude of the error we take ^ = |. It is likely that a stronger bound on 
the error is achievable, but flA.lOp is sufficient for this proof. The primary difficulty is bounding 
Cfc-i — Cfc) ill ( lA.gp . for which f lA.3|) was used crudely to obtain |Ca;-i — Cfc| < 2kMi. 

Next, (1A.4P is proved. By flA.2p . for all k > e"^ we have \yk\ > ^^(i-^). Thus by (jA.ll) . we 
can write 

yk+2 = yk-Hiyk;e)e-^^'+'^ , (A.ll) 

where ff is non- negative, bounded and has bounded derivatives. By (1A.11|) and (13. 5p . 

Vk+i = yk-{l- q[y,^,,y,])H{yk)e'^^''''^ . (A.12) 

Further applications of (lA.llI) and (13. 5p with (IA.12P lead to 

_ H{yk+2) - H{yk+i) + q[y^+^,y^]H{yk) 
1ly.,.y.,.} - J^^^ ■ (A.13) 

By Taylor expanding the various instances of H in (1A.13P about yk, we find that all 0(1) and 
0(52(1+^)) terms cancel producing, q[y^_^^,y^^2] — q[yk+2,yk\ ~ 0{e^^^). Moreover, assuming the 
vector field, (13. ip . is at least piecewise-C^, there exists a sufficiently small range of y and e values 
over which H"{yk) is bounded, which establishes f lA.4p . 

Finally, for e-independent values, |/min and i/max, ^min and /cmax are 0{e~^). Thus, by flA.4p 
and flAlQ]) with ^ = |, 

^[^^.in+-...„,J = j^^ + Oie'^) , (A.14) 
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and for any odd k with, kmin + 2 < /c < /Cmax, the difference, (i[yk,y^,_2] — 1[yk ■ +2,yk ]' ^^ ^(^0- 
Therefore, by fl3.15p . 

A 1 

9[!/n.in,?/.nax] = ^^^^g + ^(e^) , (A.15) 

which completes the proof. (Additional arguments to maximise the exponent in the error term 
of (lA.lSp are beyond the scope of this paper.) 



B Derivation of ( 15:231 ) 



li B = A^, then by stretching the negative UQ-a.xis by a factor, -j, we find that f l5.20p behaves 
similarly for negative and positive values of uq. To exploit this symmetry, let 

Z(Mo,r) = g(Mo,r) + ^g(-^,r) . (B.l) 

From fl5.20p - fl5.2ip . for uq > and r eM, Z satisfies the boundary value problem 

dz dz id^z ,^^, 

or ouq 2 ouq 

dZ 
lim Z{uo,r) = 1 , T^(0,r) = , lim Z{uo,r) = 1 . (B.3) 

r-s>-oo dUQ no-s>oo 

The unique solution to this problem is Z{uo,r) = 1. By f lB.ip . for all r, (5(0, r) = j;^Z{0,r) = 
YX4. Thus, in particular, fPno = linir-H.oo <5(0,r) = j^, as required. 



C Derivation of (15.241) 



If i3 = 1, then f l5.16p has two special properties that provide simplification. First, v{s) increases 
at a unit rate at all times. Assuming, for simplicity, v{0) = 0, we therefore have 

du{s) = \-^'' l^l \ds + dW{s), (C.l) 



and Qno = linis^.oo P[^(s) > 0] Second, we may write fIC.ip as 



du{s) = (j){u,s)ds + dW{s) , (C.2) 

where 



The key observation is that the drift in fIC.Sp is odd with respect to time, that is. 



0(n, -s) = -(/.(u, s) , (C.4) 

for all u and s. For this reason, roughly speaking, drift over negative times is cancelled by drift 
over positive times and consequently the probability that u{s) is positive approaches | as s — ?■ oo. 
Indeed (15.241) is a consequence of the following general result: 
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Lemma 3. Suppose (p is a bounded, measurable function satisfying ( C.4 ) with possibly finitely 
many values of u at which it is discontinuous. For any T > 0, for the stochastic process liC.S^) . 
F[u{T) >0 I n(-T) = 0] = i. 

The desired result, fl5.24p . follows by applying Lemma [3] to (IC.3|) and taking T — )■ 00. Below, 
Lemma [3] is proved for smooth 0. For brevity a generalisation permitting discontinuities in 
is omitted. This may be achieved by simply adding consistency conditions for the transitional 
probability density function of ( ]C.2I) at points of discontinuity. 

For f lC.2p . the probability density function of u{s), call it p{u,s), is the unique bounded 
solution to 

Here we let 

QK, So) = P [n(T) > I n(so) = ^o] . (C.6) 

By the Feynman-Kac formula [63l |59], Q is the unique, bounded, C^ solution to 

^ = -'^(-O'^o)^^ - 2^ ' ^(^o,T) = X[o,^)K) . (C.7) 

By the symmetry property, f lC.4p . the solutions to these two partial differential equations are 
related by ^(mq, Sq) = p{uo, — Sq), for all Uq and sq. Now define 

/oo 
p{u,-so)Q{u,so)du . (C.8) 

-00 

A quick manipulation shows that this function takes the value | for any value of Sq: 

/°° do 1 ("^ d 1 

-—{u,so)Q{u,so)du = - -^{Q{u,So)f du = - , (C.9) 

.00 ouo 2 J_^ du 2 

since for any fixed Sq, Q{uq, Sq) — > as Mq — > —00 and Q{uq, Sq) —)■ 1 as Mq — )• cxd. Also, by the 
strong Markov property, for any Sq, 

/oo 
p{u,SQ)Q{u,SQ)du. (CIO) 

-00 

By then setting sq = we deduce 



P [^(r) > I u{-T) = 0] = /i(0) = - , (C.ll) 



as required. 



D Derivation of (15.251) 



It easier to treat the limit ^ — t- oo. Here we show that 



lim qUAB) = \, (D.l) 
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for all B > 0. In view of (12. 6p . this is equivalent to f l5.25p . 

For fl5.16p . the limit, A — t- oo, represents infinitely strong drift for n < 0. Therefore a sample 
path from a point {u,v) with u < and t' < is immediately sent to {0,v), with probability 
1. If instead u < and v > 0, then the u- value of the solution goes to — oo, with probability 
1. Consequently, it suffices to analyse dynamics for m > and impose a reflecting boundary 
condition on the negative f-axis and an absorbing boundary condition on the positive w-axis. 
Specifically, we may assume v{s) = s, and that u{s) is governed by 



du{s) = sds + dW{s) . 
The probability density function for u{s), call it p{u, s), therefore satisfies 



dp 
ds 



dp 1 d'^p 



du 



dp 
du 



(0, s) = 2sp(0, s) , for s < 
lim p{u, s) = 6{u) , 



2(9m2 ' 

p(0, s) = , for s > 
lim p{u, s) = . 



(D.2) 

(D.3) 

(D.4) 
(D.5) 



Note, in particular, the behaviour of fl5.16p in the limit, ^ — t- oo, is independent of B. Since, 
by fl5.24p . lim^_5.oo Qnoi-A, 1) = |, we therefore have (ID.ip for all B, as required. However, to add 
to our understanding of the system for large (and small) values of A, here we derive an explicit 
expression oi p{u, s) from flD.3P - flD.5p . and use it to derive fID.ip directly. 

The problem, fID.SP - flD.Sp . with instead a refiecting boundary condition at u = for all values 
of s, was solved by Knessl in [M]. Thus for s < 0, p{u, s) is Knessl's solution: 

(D.6) 



p(u,s) = 23e-V+^«ii'(u,s) , 
where K may be expressed as an inverse Laplace transform 

I p Aif2i(n 
27ri 



Kiu.s) 



a 



Br 



Ai^ 23a 



da 



(D.7) 



In flD.7p . Ai(z) is the Airy function of the first kind, and Br denotes a suitable Bromwich contour 
in the complex plane. Throughout this exposition Bromwich contours may be assumed to be 
vertical lines at a non-negative real value. The following expression of p for s > was obtained 
by methods similar to those in [61] (for brevity the details are omitted): 



p(m. 



2 

23e" 



27ri 



,/3s 



Br 



-K{i,Q)^{u,P,Odidp 



where 



$(w,/3,0 



237rAi(23(5+/3)') 



Ai(2^/3j 
Aif2i/3 



Ai {2\{u + /?)) Bi {2\p\ + Ai {2\p\ Bi (^25 [u + /3) 
Ai {2\{i + /?)) Bi (^2^) + Ai (2M) Bi {2h{i + /3) 



(D.8) 



u <i 

u> E, 
(D.9) 
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and Bi{z) is the Airy function of the second kind. We may use flD.81) to derive (ID.ip because 



hm gno{A,B)= hm / p{u, s) du 



A^cx 



(D.IO) 



For large s, due to the nature of Brownian motion a solution to ( 1D.2I) that has not been 
absorbed at m = will very likely lie an O (52 1 distance from the value |s^ (the solution to 
f lD.2|) in the absence of the noise term). For this reason we write 



u 



5^ + As2 



(D.ii: 



and take n, s — )■ 00 with A = 0(1). 

As u — )■ 00, only the piece of (ID.9P with n > ,^ is important. From this we utilise various 
basic properties of Airy functions to obtain, for large u, 



where 



F{P) 



K(e,0)<l>(n,/3,Orfe ~ Ai(23(n + /3))F(/3) 

1 



27ri iBr Ai (2ia) Ai (2^/?) {a - /3) 



da 



subject to the restriction 



Re(a) > Re(/3) . 
Next, to ( lDl8l) with flPTTTj) and (lDn2|) . as in [M] we apply the formula, Ai(x 



valid for large \z\ bounded away from the negative real axis 

1 

1 



which yields 



1 

1 1 

27r5z5 



(D.12) 
(D.13) 
(D.14) 

e-^^ 

(D.15) 

(27ri)2 y^jBr Ai(a)Ai(6)(a - b) ^-^^ ' ^^'^^^ 

where we have substituted, a = 23a, and, b = 23/3, and hence Re(a) > Re(6). By further writing, 
a = 7 + Y + iw (for any 7 > 0), and b = — ^ + iu, we produce 

1 



p{u,s) 



A/27rs 27ri Jgr 



Therefore as s — ?■ 00, p is a scaled Gaussian and, by fID.lOp and (JD.ISP 

A^oo ZTTl 



■dadb 



00 PCX) 



lim £)no(^, -B) = "T^ / / — 7 ^ ^ 7 — ■ T^ dr duj . (D-17) 

^->cx.^°°^ ^ 47r2y_^y_^Ai(7 + f + ia;)Ai(-f + ia;)(r-i7) ^ ^ 

Since (]D.17|) is independent of the value of 7, to evaluate it we may take the limit 7—7-0. To 
achieve this we apply the following version of Plemelj-Sokhotski's formula [66]: lim^_5.o+ / —^ dx = 
mh{0) + CP. J -^ dx, (for integration over an interval containing zero), which produces 

i BUiu) °° 



1 f°° 1 
lim g^o{A,B) = -— I ^^ 



du = — 



4 Ai(icu) 



1 
2 



(D.18) 



as required. 
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